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ABSTRACT 

We present a detection of the Sunyaev-Zel’dovich (SZ) decrement associated with the Luminous 
Red Galaxy (LRG) sample of the Sloan Digital Sky Survey. The SZ data come from 148 GHz maps of 
the equatorial region made by the Atacama Cosmology Telescope (ACT). The LRG sample is divided 
by luminosity into four bins, and estimates for the central Sunyaev-Zel’dovich temperature decrement 
are calculated through a stacking process. We detect and account for a bias of the SZ signal due 
to weak radio sources. We use numerical simulations to relate the observed decrement to I 200 and 
clustering properties to relate the galaxy luminosity bins to mass. We also use a relation between 
BCG luminosity and cluster mass based on stacked gravitational lensing measurements to estimate 
the characteristic halo masses. The masses are found to be in the range ~10 13 — 10 14 /i -1 Af©, a lower 
range than has been previously probed. 
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1. INTRODUCTION 

New high resolution measurements of the sky at mil- 
limeter wavelengths are enabling the use of the ther- 
mal Sunyaev-Zel’dovich (SZ) effect as a precise cosmo- 
logical tool. The effect arises from the inverse Comp- 
ton scattering of the cosmic microwave background 
(CMB) and high energy electrons, usually occurring 
in the hot intra-cluster medium (ICM) of galaxy clus- 
ters (Zel’dovich & Sunyaev 1969; Sunyaev & ZeFdovich 
1970). For frequencies below 218 GHz, the distortion 
manifests itself as an arcminute-scale temperature decre- 
ment along the line of sight to a cluster. The amplitude 
of this decrement is nearly independent of redshift, mak- 
ing the SZ effect particularly useful for examining the 
high redshift universe. The effect can also be used to 
explore the role of baryonic physics in cluster evolution, 
as its amplitude is proportional to the thermal pressure 
of the ICM integrated along the line of sight. For a clus- 
ter in hydrostatic equilibrium, it is expected that the 
integrated SZ signal scales with total cluster mass. For a 
full description of the SZ effect, see the review articles by 
Rephaeli (1995), Birkinshaw (1999). and Caxlstrorn et al. 
( 2002 ). 

Analyses of large X-ray and optically selected clus- 
ter samples have recently illustrated the cosmological 
potential of cluster surveys (e.g., Vikhlinin et al. 2009; 
Mantz et al. 2010; Rozo et al. 2010). The emergence of 
large-area SZ cluster surveys, such as the South Pole 
Telescope (SPT; Carlstrom et al. 2009) and the Ata- 
cama Cosmology Telescope (ACT; Fowler et al. 2007; 
Swetz et al. 2010), will supplement this previous work. 
Furthermore, upcoming results from the Planck satellite 
are expected to include a large number of SZ-detected 
clusters (Bartlett et al. 2008). The SPT collaboration 
reported its first SZ detections in Staniszewski et al. 
(2009) and has since identified twenty-two cluster candi- 
dates in Vanderlinde et al. (2010). Twenty-one of these 
clusters were optically confirmed in High et al. (2010). 
The X-ray properties of the SPT sample and the de- 
rived SZ/X-ray relations are presented in a follow-up 
study (Andersson et al. 2010). The ACT collabora- 
tion presented its first cluster detections in Hincks et al. 
(2010) and has since reported on twenty-three clusters 
in Marriage et al. (2010b), all of which have been op- 
tically confirmed (Menanteau et al. 2010). The cosmo- 
logical implications of the ACT sample are discussed in 
Sehgal et al. (2010b). 

In order to fully use cluster surveys for cosmological 
purposes, independent and robust estimates of the clus- 
ter masses are needed. Several studies have explored the 
relation between the integrated SZ signal Y and total 
cluster mass M for massive clusters ( M > 10 14 Mq) (i.e., 
Benson et al. 2004; Bonamente et al. 2008; Melin et al. 
2010: Huang et al. 2010; Plagge et al. 2010). The inte- 
grated Y parameter is defined as 

Y= [ y(n)dQ, (1) 

Jn 

where y is the usual Compton -//-parameter and Q is the 
solid angle of the cluster. Assuming that thermal energy 
results solely from gravitational collapse, it is possible 
to derive self-similar scalings between observables (e.g.. 
Y) and cluster mass. Furthermore, observations show 


that the SZ-mass relation is relatively insensitive to the 
specifics of cluster physics (i.e., cooling, AGN feedback) 
(Bonamente et al. 2008; Sehgal et al. 2010b). 

In this work, we report on the stacking of a subset of 
the ACT 2009 equatorial data at the positions of Lu- 
minous Red Galaxies (LRGs) measured by the Sloan 
Digitial Sky Survey (SDSS). Various studies have exam- 
ined the clustering properties of the SDSS LRGs (i.e., 
Zheng et al. 2009; Reid & Spergel 2009) and have found 
that most of these galaxies reside in halos of typical 
mass Al ~ 10 13 — 10 14 h~ 1 AlQ. This paper describes 
the stacking process and the detection of SZ signal af- 
ter binning the sample by luminosity. Estimates for the 
characteristic halo mass of each luminosity bin are ob- 
tained using two methods: gravitational lensing mass 
measurements and an empirical model for halo bias as a 
function of luminosity (Tegmark et al. 2004; Reyes et al. 
2008; Zehavi et al. 2005) that describes the LRG sam- 
ple well, as shown in Percival et al. (2007). Finally, we 
present an initial analysis of the relation between inte- 
grated SZ signal and halo mass. Unless otherwise stated, 
we assume a flat ACDM cosmology with Q m — 0.264, F2a 
= 0.736, os = 0.80 and H a = 100/;. km Mpc -1 with 
h = 0.71 (Komatsu et al. 2009). 

This paper is organized as follows. In Section 2, we 
describe the LR.G sample and the ACT map. In Section 
3, we describe the process of filtering, and in Section 4 we 
discuss our methods for binning and stacking the ACT 
data. We present the stacking results in Section 5 and 
examine the SZ-mass scaling relation in Section 6. We 
conclude in Section 7. 

2. DATA 

In this section, we describe the subset of the ACT data 
used in this analysis (Sec. 2.1) and the LR.G sample from 
the SDSS (2.2). ' 

2.1. ACT Data 

ACT is a six-meter telescope operating in the Ata- 
cama Desert of Chile at an altitude of 5200 meters. 
The telescope site was chosen for its dry climate as 
well as the ability to observe in both the northern and 
southern skies. The telescope has three 1024-element 
arrays of transition edge sensors, one each operating 
at 148 GHz, 218 GHz, and 277 GHz. For a more 
detailed introduction to the ACT instrument, obser- 
vations, and data reduction, see Fowler et al. (2007), 
Swetz et al. (2010) and Das et al. (2010). The sci- 
ence release from the ACT 2008 survey of a region 
south of the Galactic equator includes results on the 
CMB power spectrum and related parameter constraints 
(Fowler et al. 2010; Das et al. 2010; Dunkley et al. 2010; 
Hajian et al. 2010). Results are also presented on 
compact millimeter sources (Marriage et al. 2010a) and 
clusters (Marriage et al. 2010b; Menanteau et al. 2010; 
Sehgal et al. 2010b). The present study uses data at 
148 GHz from a subregion of the 2009 equatorial sur- 
vey. The subregion covers 268 square-degrees and lies 
between right ascensions 21 h 20 m arid 03 h 40 m and decli- 
nations — ORSO' and 01°30'. This region is coincident 
with the SDSS Stripe 82, which contains a rich multi- 
wavelength data set. Figure 1 is a map of the sensitivity 
across the subregion of the study, along with the loca- 
tions of the LRGs. The map has been match-filtered for 
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FlG. 1. — Sensitivity map with LRG locations. The map shows the sensitivity over the subset of the ACT 2009 148 GHz equatorial data 
considered for this study. The gray-scale encodes the noise rms in pK of a map match- filtered for cluster detection. The median noise in 
the map is 34 pK in CMB temperature units. 


cluster detection as described in Section 3. The median 
rms in the filtered map is 34 uK in CMB temperature 
units. 

2.2. LRG Sample 

The LRG sample used in this study is from the SDSS. 
a survey that lias performed five-band ( ugriz ) pho- 
tometry (Fukugita et al. 1996) using a specialized wide- 
field camera (Gunn et al. 1998) and multi-object spec- 
troscopy using a pair of spectrographs. The survey 
has imaged one quarter of the sky at high Galactic 
latitude and conducted spectroscopic follow up obser- 
vations of approximately one million of the detected 
objects (Eisenstein et ai. 2001; Richards et al. 2002; 
Strauss et al. 2002). All data are processed automati- 
cally by pipelines that detect and measure photometric 
properties of sources (Lupton et al. 2001). The SDSS 
has had seven major data releases (Abazajian et al. 2003, 
2004. 2005; Adelman- McCarthy et al. 2006, 2007, 2008; 
Abazajian et al. 2009). 

The LRG catalog used in the current work is drawn 
from the sample prepared by Kazin et al. (2010), which 
is taken from the seventh data release (DR7) of the SDSS. 
The data set was first, developed by Eisenstein et al. 
(2001) and is publicly available. 1 The DR7 LRG sam- 
ple contains ~ 110,000 galaxies and extends to redshift 
2 ~ 0.5. The subset used in this study contains 2681 
LRGs that fall within the 268 square-degree ACT map. 
The subset has a spectroscopic redshift range 0.16 < 2 < 
0.47 and has rest frame ("/-band absolute magnitudes 
—23.2 < M g < —21.2. K corrections (Blanton & Roweis 
2007) have been applied to the galaxies. The SDSS cat- 
alog of LRGs serves as a good tracer of large scale struc- 
ture, as LRGs probe a larger volume than other tracers 
since they are more luminous than most galaxies. They 
also exhibit a distinct 4000 A break in their spectral en- 
ergy distributions, making redshift determination easier 
than for other galaxies. 

3. FILTERING THE ACT DATA 

To maximize the SZ signal of the clusters asso- 
ciated with the LRGs. we employ a matched fil- 

1 http://cosmo.nyu.edu/~eak306/SDSS-LRG.html 


ter (Haehnelt & Tegmark 1996; Herranz et al. 2002; 
Melin et al. 2006). The filtering scheme used here closely 
follows Marriage et al. (2010a), with exceptions rioted 
below. At a position x in the map, we model the temper- 
ature measured as the sum of the SZ signal plus the other 
components in the map (CMB fluctuations, atmospheric 
fluctuations, instrument noise and point sources): 

5T(x) = J2 STo,ibi(x - Xi) + h T'other)#), (2) 

i 

where ST 0t i is the peak amplitude of the i th LRG and 
bi is the unit-normalized ACT 148 GHz beam function, 
taken to be isotropic with a FWHM of 1.4' (Hincks et al. 
2010). We first mask the brightest sources in the map, 
which amount to 146 objects. Only pixels associated 
with a source with at least SNR > 5 are masked. These 
are predominantly radio sources (Marriage et al. 2010a). 
None of these sources lies within 0.5' of an LRG po- 
sition. Before applying the matched filter, we multi- 
ply the map, pixel-wise, by the number of observations 
per pixel normalized by the m aximum number of obser- 
vations in any single pixel, y/N Qbs (x)/N obs , max . This 
weighting accounts for local changes in the amplitude of 
the white noise and is equivalent to down-weighting pix- 
els with a high white noise rms. This weight scheme is 
also used during the stacking process when determining 
the matched filter decrement 6T 0 (see Section 4.2). Next, 
the map is filtered in Fourier space using a matched filter: 

= JWjs^ e ^s^ L 

Jb*(k')\ST other (k')\- 2 b(k')dk' 

Since the power spectra of the atmosphere and the CMB 
scale as l ~ 4 for l > 1000. the filter is approximately a 
fourth derivative operator convolved with the beam. The 
choice of the beam function b for the spatial profile of the 
cluster differs from the method used in Marriage et al. 
(2010b), which employs a /3-model for the cluster pro- 
file. By comparisons with simulations, we expect that 
the ACT 148 GHz beam function will be a good match 
for the spatial profile of clusters and hence adopt it in 
the matched filter. It is particularly effective at reducing 


4 


Hand efc al. 


TABLE 1 

Luminosity Bins and Halo Masses for Stacked LRGs 


Bin 

AW 

(Lo.lr) 

10 lo i G 

£■ 0.1 r Range 
10 10 L(o 

%> 

A'-boo p & 
10 I4 1Rq 

AAooph 

10 14 A/ o 

M 2 oop c 

1 

26 

30.1 

25.9 - 44.9 

0.39 

2.74 ± 0.74 

5.20 ± 1.23 

2.59 ± 0.87 

2 

60 

23.6 

21.8 - 25.3 

0.40 

1.78 ± 0.13 

3.57 ± 0.23 

2.05 ± 0.43 

3 

674 

16.7 

13.8 - 21.8 

0.37 

1.00 ± 0.22 

2.17 ± 0.41 

1.43 ± 0.31 

4 

1921 

10.0 

1.1 - 13.8 

0.30 

0.43 ± 0.14 

1.07 ± 0.29 

0.82 + 0.22 


“Determined from halo bias using 6* = 0.99. 

^Determined from halo bias using 6* = 1.30. 

“Determined from lensing measurements in Reyes et al. (2008). 

scatter from CMB noise, to which larger filters are more 
sensitive. 

Rather than explicitly modeling the noise contribution 
8T 0 ther(*)i tve use the original input map with sources 
masked for the power spectrum of ST otl i iev (k) in equation 
3. As the model for 5T a ther(*) is not well-conditioned, 
we smooth !<5T ot her(fc)l~ 2 with a Gaussian filter (a = 
5 Fourier space, pixels). We also set outlying Fourier 
space pixels with amplitude greater than 10 times the 
median pixel value to the median value. We also in- 
clude an additional low-t filter that grows from zero at 
i = 0 to unity at £ = 1200 as sin 5 (7r£ /2400) . This filter 
accounts for any modification of CMB statistics due to 
multiplying by the non-uniform weight map and makes 
our result more robust to any potential non-Gaussianity 
in the large-scale noise that might arise from residual at- 
mospheric 1 noise in the maps. After finishing the filtering 
process, we do not include any sources that fall within 
10' of the map edge in order to mitigate noise from edge 
effects. This reduces the map size from 284 square de- 
grees to 268 square-degrees. The entirety of the LRG 
sample discussed above (2681 objects) falls within the 
bounds of the reduced map, leading to roughly 10 LRGs 
per square-degree. 

4. STACKING METHODS 

While ACT can easily detect massive clusters (M > 
10 15 Mq) (Marriage et al. 2010b), the CMB fluctuations 
and detector noise limit the direct detection of less mas- 
sive clusters. Since the mean of both the CMB fluc- 
tuations and the atmospheric noise should be zero, we 
stack the LRG cluster submaps to enable the detection 
of weaker SZ signals. 

4.1. Luminosity Bins 

Before stacking the LRGs, we bin the sample by r- 
band luminosity K corrected to z = 0.10, Lo.ir- Table 1 
shows the four luminosity bins used in the analysis. This 
table also lists mass estimates for each bin, as described 
in Section 5.2. The luminosities have been computed us- 
ing the Galactic extinction-corrected (using maps from 
Schlegel et al. 1998) Petrosian magnitudes from the DR 7 
CAS and their measured errors. They have been K cor- 
rected to z = 0.10 using the kcorrect.v4_2 software 
package (Blanton & Roweis 2007). There have been no 
added evolutionary corrections. We also refer to the K 
corrected a 1 r-band absolute magnitude of an LRG as 
Following Blanton et al. (2003), 

Mo.ir ~ m r — DM (z, fi m , fL\, h) — Ko.i rr (z), ( 4 ) 


where DM(z,Sl tn ,£l/{,h) is the distance modulus in a 
(fi m = 0.3, Da = 0.7 , h = 0.71) cosmology and K 0Arr (z) 
is the K correction from the r-band of a galaxy at 
redshift 2 to the 0ll r-band. Luminosities are reported 
in solar luminosities, where we have used M 0l o.i r = 
4.76 (Blanton et al. 2003). These definitions of luminos- 
ity and magnitude will facilitate the estimation of halo 
masses for our binned data in Section 5.2. 

4.2. Algorithm 

We stack filtered submaps with area 10' x 10' centered 
on each LRG in a given luminosity bin. and then esti- 
mate the matched filter decrement ST 0 from the central 
temperature value in the stacked map. Recall that the fil- 
tered map is weighted by the number of observations per 
pixel normalized by the maximum number of observa- 
tions in any single pixel. To minimize the errors induced 
by the 0.5' pixel size of the ACT maps (Marriage et al. 
2010a), the pixel size of the submap centered on each 
LRG is decreased using Fourier interpolation. For the 
present analysis, the pixel size is decreased to 0.03125' 
and the submaps are recentered. Before determining 
5T 0 , the stacked map for a given bin is convolved with 
the ACT 148 GHz beam function b to recover decre- 
ment information from the surrounding pixels. Lastly, 
wc sum over the five pixels within 4" of the centroid of 
the stacked map and divide by the sum of the weights to 
compute 5T 0 . We compute the weight for each submap 
in the stac k before decrea sing the pixel size by using the 
value of v/Aobs/jVobs.max at its central 0.5' pixel. 

5. STACKING RESULTS 
5.1. Matched Filter Decrement 8T 0 

Tables 2 and 3 contain estimates for the matched filter 
decrement ST a obtained for each of the four luminosity 
bins. Due to radio contamination (see discussion below), 
we report results for the full LRG sample as well as a 
radio-quiet, subset, which excludes approximately 10% of 
the LRG sample. The value for 5T a reported in Tables 
2 and 3 is the temperature decrement in the ACT map, 
recovered from the matched filter. Figure 2 shows 7' x 7' 
stacked submaps for each luminosity bin for both LRG 
samples. Stacked maps obtained when using random 
LRG positions are also shown, and it is apparent that 
there is no clear detection in these maps, as expected. 

In addition to the central peak in the stacked images, 
several other significant peaks are evident, particularly 
the deep cold spot at the upper edge of the stacked map 
for Bin 2. which has a depth of around 24 pK and is 
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TABLE 2 

Central SZ Temperature Decrement and Integrated SZ Signal for All Stacked LRGs. Binned by Lo.ir 


Bin 

■^bin 

ST 0 

pK 

SNR 

*200p 

aremin 2 

V200P <T a E(z) 273 
Mpc 2 

1 

26 

-27.5 ±11.3 

2.4 

(1.5 ± 0.5) xlO" 4 

(1.1 ± 0.4) xlO-' 5 

2 

60 

-9.9 ± 5.3 

1.9 

(4.7 ±2.1) x 10~ 5 

(3.4 ± 1.8) x 10 -6 

3 

674 

-4.7 ± 1.5 

3.1 

(1.8 ± 0.5) xlO- 5 

(1.3 ± 0.4) xlO -6 

4 

1921 

-1.6 ± 1.0 

1.6 

(3.7 ± 4.6) xl0~ B 

(3.8 ± 2.9) xlO" 7 


TABLE 3 

Central SZ Temperature Decrement and Integrated SZ Signal for Stacked Radio-quiet LRGs. Binned by Lo.ir 


Bin 

N'bin 

ST 0 

mK 

SNR 

^200/5 

aremin 2 

Y 20 op rf 2 E(z )-*/ 3 
Mpc 2 

1 

21 

-28.2 ± 12.9 

2.2 

(1.5 ± 0.6) xlO -4 

(1.1 ± 0.5) xlO" 5 

2 

51 

-10.6 ± 4.9 

2.2 

(4.2 ± 1.9) xlO -5 

(3.7 ± 1.7) x 10 -B 

3 

587 

-6.1 ±1.6 

.3.8 

(2.4 ± 0.5) xlO” 5 

(1.7 ± 0.5) X10-® 

4 

1732 

-2.1 ± 1.0 

2.1 

(7.8 ± 6.9) x 10-® 

(6.0 ± 4.5) x 10“ 7 


Note. — The estimates for Y 2 oop are determined from 5T a using microwave sky simulations (Sehgal et al. 2010a), They correspond 
to the halo masses M 2 oop derived from considering luminosity-dependent halo bias with 6* — 0.99 (Table 1). The uncertainty in STo is 
determined from bootstrap resampling. The central decrements and associated errors are measured in CMB temperature units. 


significantly deeper than the stacked cluster at the cen- 
ter. For this stacked map. the pixel noise is around 5 
|xK, making the upper-edge cold spot a 5a discrepancy. 
Examination of the 60 individual maps in the stack re- 
veals that two maps contain decrements of < 150 pK at 
the location of the cold spot. In addition, 11 additional 
map sections show decrements larger than 50 uK at the 
approximate position of the upper-edge cold spot. This 
large number of decrements in similar positions relative 
to the LRG’s in our catalog clearly is not, a random phe- 
nomenon, but rather reflects the fact that galaxy clusters 
and groups are strongly clustered on the sky. 

The errors on 5T 0 presented in Tables 2 and 3 are 
computed for each bin using using bootstrap resampling. 
These errors measure the scatter in the properties of in- 
dividual LRGs. Figure 3 shows the distribution of the 
mean of ST 0 for each of the four luminosity bins for the 
full and radio-quiet LRG samples. The distribution of 
each bin is well-approximated by a Gaussian distribu- 
tion. 

We use numerical simulations (Sehgal et al. 2010a) to 
relate the matched filter decrement to the integrated Y 
parameter (equation 1). In an unfiltered map, where the 
cluster is larger than the beam, Y ~ SToO^., where 0 C 
is the angular size of the cluster. For a cluster of fixed 
physical size, Y oc ST 0 d ^ where d A is the angular diam- 
eter distance to the cluster. Because the filter acts as a 
V 4 operator and removes large-scale signals in the map, 
the expected relationship for nearby clusters of fixed size 
is likely even steeper. However, rather than rely on this 
approximate scaling to weight the observations, we apply 
the same filter used on the ACT data to the simulated 
maps and compute 5T 0}i for each of the clusters. The 
simulated maps are first convolved with the ACT 148 
GHz beam function. We then fit a simple model to the 
Y — ST relation based on the angular diameter distance 
scaling: 

Yesui = C ST. j z'~ fi . (5) 

where Zj is the redshift and Y est j serves as an estimate 
for the SZ signal within a disk of radius i? 20 o for the 2 th 


cluster. We define R 200 as the radius within which the 
average density is equal to 200 times the mean matter 
density of the universe, p(z ) = p CT it(z — 0)f2 m (l + z) 3 . 
The relationship between Y and ST a is obtained for the 
mass range of each luminosity bin. Here the mass range 
is represented by the la uncertainty in the central AT oop 
value, listed in Table 1. Due to systematic uncertainties 
in determining halo masses, we perform this analysis for 
three separate groups of mass bins (see Section 5.2). 

The best fit values and corresponding errors for C and 
P are summarized in Table 4 . The uncertainties reported 
are the formal la errors computed using Levenberg- 
Marquardt least squares minimization. The correlation 
between Y es t and Spoop is shown in Figure 4. Figure 5 
shows the redshift dependence inherent in the relation 
between Y and ST 0 , which motivated our model for Y e$t 
in equation 5. Figures 4 and 5 show the results for the 

TABLE 4 

Best-Fit Parameters for the V200P - ST 0 relation, obtained 
FROM SIMULATION AND PITTED USING EQUATION 5 


^Vbin 

Mass Range 
10 I4 AT G 

c 

aremin 2 /K 

0 

Halo Bias using 6* = 0.99 

1103 

2.00 - 3.48 

1.37 ± 0.05 

1.25 ± 0.03 

1245 

1.65 - 1.91 

1.30 ± 0.05 

1.18 ± 0.03 

10321 

0.78 - 1.22 

1.08 ± 0.02 

1.16 ± 0.02 

15116 

0.29 - 0.57 

0.92 ± 0.05 

1.1-1 ± 0.04 

Halo Bias using 6* = 1.30 

r 382 

3.97 - 6.43 

1.47 ± 0.11 

1.53 ± 0.07 

291 

3.34 - 3.80 

1.42 ± 0.11 

1.37 ± 0.07 

2152 

1.76 - 2.58 

1.35 ± 0.04 

1.24 ± 0.02 

11274 

0.78- 1.36 

1.09 ± 0.02 

1.19 ± 0.02 

Weak Lensing Measurements 

3127 

1.72 - 3.46 

1.37 ± 0.03 

1.30 ± 0.02 

590 

1.62 - 2.48 

1.21 ± 0.07 

1.24 ± 0.05 

5408 

1.12 - 1.74 

1.16 ± 0.03 

1.20 ± 0.02 

11274 

0.60 - 1.02 

1.12 ± 0.03 

1.08 ± 0.02 


Note. — The three groups of mass ranges correspond to the 
three separate methods that are used to obtain halo masses for the 
LRG sample. 
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mass bins corresponding to 6* — 0.99 in Table 1 . 

The values for Ymop are found for the j th luminosity 
bin by summing over each LRG in the bin as follows. 


TABLE 5 

Matched Filter Decrement ST 0 for Various Systematic 
Tests of the ACT Data 




200 p,j 


V; V ' Wi Cj 5T 0 j z, 01 


E 


Nj 

Wi 


( 6 ) 


where <5T 0> » is the matched filter decrement recovered 
from the filtered ACT map and Nj is the number of 
LRGs in the j th luminosity bin. The weight uy is re- 
lated to the number of observations per pixel (see Sec- 
tion 4.2). Values for 1200,5 are reported in Tables 2 and 
3. These values correspond to the halo masses Ahoop 
derived from considering luminosity-dependent halo bias 
with 6* = 0.99. The error in represents the uncer- 
tainty after propagating errors in ST 0 , /l and C . These 
formal errors do not include the systematic uncertain- 
ties in our model which assumes that the gas profiles in 
Sehgal et al. (2010a) are accurate representations of the 
profiles for typical clusters. 

Table 5 contains the matched filter decrements 5T 0 ob- 
tained when stacking various subsets of the LRG sample. 
These subsets are designed to test various systematic as- 
pects of the ACT data. The first four tests in Table 5 
probe spatial discrepancies in the data by splitting the 
ACT map by right ascension and declination. We con- 
clude that the signal is consistent throughout the map. 
We also examine variations in white noise amplitude by 
splitting the sample spatially by the number of observa- 
tions JVohs per pixel. We determine the median value of 
A’ 0 bs per pixel and then stack subsets of LRGs that lie 
above and below this median value. We conclude that 
the signal is present for both high and low noise bins in 
the map. 

Radio source contamination at 148 GHz appears to 
have a significant effect. Using the SDSS DR7 CAS, we 
identified LRGs with radio counterparts within 2" in the 
1.4 GHz VLA FIRST Survey (Becker et al. 1995) and 
define these LRGs to be radio-loud. There are 290 radio- 
loud LRGs in the sample: 5/26 in Bin 1, 9/60 in Bin 2, 
87/674 in Bin 3 and 189/1921 in Bin 4. Contamination 
by FIRST radio sources is potentially a concern for ~ 
10% of the LRG sample. We also compute the matched 
filter decrement for each bin excluding the 290 radio- 
loud LRGs. While the stacked radio-quiet LRGs show a 
significant (> 4cr) decrement, the radio-loud LRGs show 
clear evidence of radio source contamination. For the 
rest of this analysis, we report results for both the full 
and radio-quiet samples (see Tables 2 and 3). 


5.2. Estimates for Halo Mass M 200 p 

We take two different approaches to estimating the 
characteristic mass associated with halos of a given lu- 
minosity. The first approach is to use measurements of 
the halo bias, defined as the relation between the galaxy 
power spectrum P gg (k) and the linear matter power spec- 
trum Plin(k), 


Pgg(k) = b(k) 2 Pu n (k). (7) 

On the large scales probed by the LRGs, the bias b(k) 
is expected to be nearly scale-independent, but strongly 
dependent on luminosity. Tegmark et al. (2004) and 


Test 

-^bin 

STo 

(mK) 

Top 

1295 

-2.51 

± 

1.16 

Bottom 

1386 

-3.13 

± 

1.10 

Left 

1406 

-2.94 

± 

0.99 

Right 

1275 

-2.67 

± 

1.35 

High N obs 

1455 

-2.15 

± 

0.98 

Low N a bs 

1226 

-3.85 

± 

1.36 

Radio- loud 

290 

3.27 

± 

2.65 

Radio-quiet 

2391 

-3.56 

± 

0.84 


NOTE. — The central decrements and associated errors are mea- 
sured in CMB temperature units. 


Zehavi et al. (2005) have fitted an empirical model for 
galaxy bias relative to the bias of L* galaxies 1 : 

- = 0.85 + 0.15— + 0.04(M* - M 0 . m), (8) 

0* L/* 

where Af* — 51ogxo(/t) = —20.44 (Blanton et al. 2003). 
For h — 0.71, M* = —21.18. Both M* and L, are in the 
01 r-band. Percival et al. (2007) have shown this model 
for relative bias to be a reasonable fit for the LRG sample 
considered in this paper. 

There exists some uncertainty in the bias normaliza- 
tion factor, 6*. Because this uncertainty can lead to large 
differences in halo mass, we will report mass estimates for 
two estimates of 6*. The values used here are 6* = 0.99 
(Seljak et al. 2005) and 6* = 1.30 (Reid et al. 2010). 

Tinker et al. (2010) use numerical simulations to de- 
rive a relationship between halo mass and bias. They 
identify halos in their cosmological simulation using the 
spherical overdensity algorithm to define the halo mass: 

M a = ~ttR 3 a p(z)A, (9) 

where p(z) is the mean matter density of the universe at 
redshift z. This analysis will use A = 200 and denote 
the halo mass by iVfy oop- Tinker et al. (2010) estimate 
the relation between halo bias and M2oop- 

b(y) = 1 - A s + 0.183J/ 1 - 5 + Cv 2A , (10) 

v° + dy 

where v represents the ‘‘peak height” of the density field 
and is given by v = Sc/v^'hoop), S c is the critical density 
for collapse, and a is the linear matter variance on the 
scale of the halo, i.e., R = (SM/Anp) 1 ^ 3 . We use z = 0 
to determine the matter variance a in order to compare 
b(v ) to b/b*. In all calculations, we use 8 C = 1.686. The 
parameters in equation 10 are: 

A— 1.0 + 0.24y exp[— (4/t/) 4 ], 

B = 0.44y -0.88, 

C = 0.019 + 0.1072/ + 0.19 exp[-(4/t/) 4 ], 
where y = log 10 (200). 

1 Lt, is defined with respect to the Schechter (1976) luminosity 
function for SDSS galaxies (see Blanton et al. 2003) 
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Using equations 8 and 10, we obtain estimates for halo 
mass for each luminosity bin. These values are given in 
Table 1. We determine mass from luminosity for each 
LRG in a given bin, and report M 200 p as the mean of 
these values. The error on M^oop is estimated by the 
standard deviation of this mean, within each bin. 

Stacked gravitational lensing measurements are an al- 
ternative method of estimating the halo mass-galaxy lu- 
minosity relation. Reyes et. al. (2008) stack a sample of 
clusters from the SDSS maxBCG catalog (Koester et al. 
2007a, b) and use weak gravitational lensing measure- 
ments to investigate the relation between cluster mass 
and various optical tracers. Since LRGs trace massive 
halos (Ho et al. 2009), the maxBCG catalog should rep- 
resent. a population similar to the LRG sample used in 
this work. Noting that Reyes et al. (2008) bin their clus- 
ter sample by °' 25 r-band BCG luminosity, we compute 
the 0,25 r-band luminosity for each LRG and determine 
halo masses using equation 15(c) of Reyes et al. (2008). 
Table 1 lists these mass values, as well as those obtained 
from analysis of halo bias. The range in these three 
masses represents some of the systematic uncertainties 
in relating LRG luminosity to halo mass. 

6. SZ-MASS SCALING RELATION 

In this section, we present an analysis of the correlation 
between cluster mass and integrated SZ signal. Specif- 
ically. we compare Tjoop to the cluster mass Mzaap for 
the luminosity/mass bins in Table 1. For self-similar evo- 
lution, analytical theory (e.g., Komatsu & Seljak 2001) 
predicts a simple relationship between the integrated Y 
parameter and cluster mass. Assuming virial and hydro- 
static equilibrium, the cluster gas temperature can be 
related to M and E(z) through (e.g.. Bryan & Norman 
1998) 

T ozM 2/3 E(z) 2/3 , (11) 

where for a flat ACDM cosmology 

E(z) = [n m (i + zf + n A } l /\ (12) 

Under assumptions of an isothermal ICM, 

^200p oc d A 2 {z) Mooop T, (13) 

where d A (z) is the angular diameter distance to the clus- 
ter. Combining the above equations leads to the expec- 
tation for self-similar evolution 

>200, X d~ 2 M^oop E{zfl\ (14) 

Anticipating this scaling relation, we have plotted 
Y 2m d? A E- 2li against the halo masses Afaoop in Figure 
6. We show the three separate halo mass estimates for 
a given bin to illustrate the uncertainties that exist in 
converting LRG luminosity to halo mass. The results for 
both the full (Table 2) and radio-quiet (Table 3) LRG 
samples are shown in this figure. The expected Y — M re- 
lation as determined from numerical simulations is shown 
as a solid line and is given by: 

Y 2 oopd 2 A E- 2 ' 3 ‘= Iff* (yflfg:) • (15) 


where the best fit parameters are a = 1.76 and 7 
= -5.74. These parameters are obtained from the 
microwave sky simulations described in Sehgal et al. 
(2010a), where halos of mass M200P > 2.82 x 10 13 M 0 
have been considered. Furthermore, the redshift range 
of the clusters is restricted to 0.15 < z < 0.50, and a 
cylinder definition of Y is used. Tables 2 and 3 give 
the values of Y2oopd A E~ 2 / 3 and associated errors for 
the full and radio-quiet samples. These estimates are 
determined using a summation similar to equation 6, 
where the redshift of each LRG is used to calculate 
d A {z) and E(z). They correspond to the masses derived 
from considering luminosity-dependent halo bias with 6* 
= 0.99. 

7. DISCUSSION 

By stacking LRGs, we have obtained estimates of the 
SZ signal as function of mass for groups and clusters. The 
stacking analysis enabled the detection of an SZ signal 
for clusters with masses as low as ~3 x 10 13 A/q, and the 
amplitude of the signal is consistent with our theoretical 
estimates. Using masses derived from halo bias, the slope 
of the inferred Y—M relation is consistent with expected 
results. However, when using stacked weak gravitational 
lensing measurements to estimate the halo masses, the 
slope appears to be slightly steeper than the predicted 
value. Agreement with the overall normalization of the 
Y — M relation depends upon our method of estimating 
the halo mass from luminosity, which is subject to large 
uncertainties. 

We have detected radio source contamination of the 
SZ signal at 148 GHz for ~10% of the LRGs. This sug- 
gests the possibility that infrared sources may also have 
a contamination effect on our measurement of the SZ sig- 
nal for the LRG sample. Mittal et al. (2009) found that 
the likelihood for a cluster to host a radio-loud BCG 
is a strong function of the central cooling time, in the 
sense that all strong cool-core clusters contain a radio- 
loud BCG. Thus, our radio-quiet LRG sample may select 
against clusters/ groups with a strong cool core. However, 
in his analysis of a large sample of clusters and groups at 
0 < 0.1, Sun (2009) did not find any radio-loud central 
galaxies in groups with strong cool cores. If this is true 
at higher redshift, our radio-quiet LRG sample would be 
representative of all groups and clusters. 

Another significant source of systematic error is the 
conversion of the measured decrement 5T 0 to the inte- 
grated SZ signal Y. Our approach relies on simulations 
from Sehgal et al. (2010a) to estimate the effect of the 
matched filter on smoothing the cluster signal. Future 
work will explore this uncertainty by using simulations 
with different models of duster gas physics. 

This work is a first application of a novel stacking 
method for searching for the SZ effect in clusters by ap- 
plying a full optical catalog to a CMB survey. Both data 
sets are rich and have multiple components, with pos- 
sible correlations between them, that we will be sorting 
out as the data and our understanding improves. Future 
ACT analysis of the SZ effect from galaxy clusters will 
include data from multiple seasons and millimeter bands 
and will include the new SDSS DR8 analysis. Our work 
shows that the presence of radio sources coincident with 
LRGs affects the 148 GHz SZ signal. Our maps show 
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hints that neighboring clusters are also affecting the sig- 
nal. We plan to explore these effects in subsequent stud- 
ies. 
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(a) Full LRG Sample 


(b) Radio-quiet LRG Sample 


(c) Random Sample 


Fig. 2. The stacked filtered temperature maps for each luminosity bin for the (a) full LRG sample, (b) radio-quiet LRG sample and 
(c) a randomly selected sample. Each map is 7' x 7' and is in CMB temperature units (pK). The values for ST C and the number of maps 
in each stack are given in Tables 2 and 3 for (a) and (b), respectively. The radio-quiet sample excludes the 10% of the full LRG sample 
that has been defined as radio-loud (see Section 5.1). The stacked random maps serve as a null test of the data, as the SZ signal should 
be consistent with zero for these maps. The number of maps in each stack for the random sample is equal to that of the full sample. The 
circles in each submap mark the characteristic size of a cluster for a given luminosity/ mass bin, in terms of #200- The values for #200 have 
been determined from simulation. The crosses mark the center of each map. 
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simulated Y 200 - (10 3 arcmin 2 ) simulated Y 200 ^ (10 3 arcmin 2 ) 


Fig. 4. — The correlation between ST 0 and for each luminosity/mass bin, as determined from simulations described in Sehgal et ai. 

(2010a). The simulation values for >200p of the i th cluster are fit to the model Y es t } i — C ST 0> i z~^ . The four subplots correspond to the 
luminosity bins used in the LRG stacking analysis. The mass range for each bin corresponds to the mass estimates determined from using 
halo bias, with 6* —0.99 (Table 1). 
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simulated Y 200 p (10 3 arcmin 2 ) simulated Y 2(IOp ( 10 3 arcmin 2 ) 



simulated Y 200i5 (10 3 arcmin 2 ) simulated Y 200 ^ (10 3 arcmin 2 ) 


Fig. 5, — The redshift dependence of the relation between 5To and Y 2 oop for each luminosity /mass bin, as determined from simulations 
described in Sehgal et al. (2010a). The redshift ranges are represented as follows: 2 < 0.26 (blue squares), 0.26 < z < 0.37 (red circles), 
and 2 > 0.37 (green triangles). The four subplots correspond to the luminosity bins used in the LRG stacking analysis. The mass range 
for each bin corresponds to the mass estimates determined from using halo bias, with 6* = 0.99 (Table 1). 
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FlG. 6. — SZ -- Mass Correlation. Yzoop estimates are plotted against cluster masses for each luminosity bin for both the full (Table 2) and 
radio-quiet (Table 3) LRG samples . The three sets of points shown use three separate mass estimates for the four LRG luminosity bins. 
The black circles correspond to mass estimates derived from analysis of halo bias using 6* = 0.99 (Seljak et al. 2005). The red triangles 
correspond to mass estimates derived from analysis of halo bias using 6* = 1.30 (Reid et al. 2010). Finally, the blue squares correspond to 
mass estimates derived from weak gravitational lensing measurements (Reyes et al. 2008). The range in these masses represents some of 
the systematic uncertainties in converting LRG luminosity to halo mass. The solid line shows the expected model for the Y — M relation 
(equation 15), as determined from the microwave sky simulations of Sehgal et al. (2010a). 





